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In addition to a smooth component of WIMP dark matter in galaxies, there may 
be streams of material; the effects of WIMP streams on direct detection experiments 
is examined in this paper. The contribution to the count rate due to the stream cuts 
off at some characteristic energy. Near this cutoff energy, the stream contribution 
to the annual modulation of recoils in the detector is comparable to that of the 
thermalized halo, even if the stream represents only a small portion (~5% or less) 
of the local halo density. Consequently the total modulation may be quite different 
than would be expected for the standard halo model alone: it may not be cosine-like 
and can peak at a different date than expected. The effects of speed, direction, 
density, and velocity dispersion of a stream on the modulation are examined. We 
describe how the observation of a modulation can be used to determine these stream 
parameters. Alternatively, the presence of a dropoff in the recoil spectrum can 
be used to determine the WIMP mass if the stream speed is known. The annual 
modulation of the cutoff energy together with the annual modulation of the overall 
signal provide a "smoking gun" for WIMP detection. 



I. INTRODUCTION 



The Milky Way, along with other galaxies, is well known to be encompassed in a massive 
dark matter halo of unknown composition. Leading candidates for this dark matter are 
Weakly Interacting Massive Particles (WIMPs), a generic class of particles that includes 
the lightest supersymmetric particle. Numerous collaborations worldwide have been search- 
ing for these particles. Direct detection experiments attempt to observe the nuclear recoil 
caused by these dark matterparticles interacting with nuclei in the detectors. These experi- 
ments include DAMA/Nal Q, DAMA /LIBRA ||, NAIAD 0, CDMS H, EDELWEISsB , 
ZEPLIN EL XENON 0, DRIFT CRESST 0113, SIMPLE HI, PICASSO Q, 

COUPP |l4j . and many others. It is well known that the count rate in WIMP direct detec- 
tion experiments will experience an annual modulation pj p| as a result of the motion of 
the Earth around the Sun: the relative velocity of the detector with respect to the WIMPs 
depends on the time of year. 

The actual signals in a dark matter detector, including the modulation, depend on the 
distribution of WIMPs in the galaxy. It is commonly believed that a majority of the WIMPs 
have become thermalized into a smooth halo distribution. In the simplest model, the halo 
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is a spherically symmetric, non-rotating isothermal sphere of WIMPs |16[. The annual 
modulation of non-standard halos, such as anisotropic models, is discussed in Refs. [l7, 18. 

However, galaxy formation is a continual process, with new material still being accreted 
through, e.g., absorption of dwarf galaxies such as the Sagittarius dwarf galaxy (22l. I23L E3| . 
The result is that the halo can contain a non-trivial amount of substructure such as clumps 
and tidal streams of material. The presence of such substructure is supported by iV-body 
simulations of galaxy formation (see, e.g., Refs. jjjl \26]). Unlike the virialized component 
of the halo, a clump or stream of material would result in a "cold" flow of WIMPs through 
a detector: the velocity dispersion is small relative to the typical speed with respect to the 
Earth, so that the WIMPs are incident from nearly the same direction and with nearly the 
same speed. Alternative models of halo formation, such as the late-infall model [27|, |2Jj |29j 
(more recently examined by Sikivie and others 0, EH El, EH EH E^), also predict cold 
flows of dark matter. Any such streaming of WIMPs (we will henceforth use "stream" 
to imply any cold flow) will yield a significantly different modulation effect than that due 
to a smooth halo. A stream actually provides two types of modulation: a modulation in 
the overall signal and a modulation in some cutoff energy above which counts due to the 
stream are not observed. Together, these two types of modulation can yield a "smoking 
gun" for WIMPs. The modulation in the presence of one or more streams is discussed in 

Refs. EES El, EI S3- 



Here, we examine how various parameters describing a stream affect the modulation 
signal; the parameters examined are the stream direction, speed, density, and dispersion. We 
expect the stream density (and hence contribution to the recoil rate) to be small, C(few%) 
that of the isothermal Halo, and yet find that the stream can have significant effects on the 
annual modulation. We note the following interesting modifications to the overall annual 
modulation when we add the effects of the stream into the isothermal halo: 

• The combined modulation generally does not have a cosinuoidal variation with time, 
even if the modulation of each individual component does; 

• The combined modulation need not be time-symmetric, even if each individual com- 
ponents is; and 

• The minimum and maximum recoil rates need not occur 0.5 years apart. 

• Near the cutoff energy of the stream (above which it no longer contributes), the stream 
contribution to the annual modulation is comparable to that of the halo, even for a 
2% stream density (relative to the background halo). The drastic effects of the stream 
near the cutoff energy can be seen in Figure ^ (for the modulation shown at 39 keV). 

• The annual modulation of the cutoff energy together with the annual modulation of 
the overall signal provide a "smoking gun" for WIMP detection. 

As a consequence of the second-to- last point, it is likely that the existence of a stream will 
be identified near its cutoff energy. Since a stream's effects are mild except near E c (and 
essentially non-existent well above E c ), the presence of a stream should not interfere with 
using the modulation to describe the background distribution, SHM or otherwise. 

In Section |Hj we will begin by reviewing direct detection techniques for WIMPs and 
the time-dependent WIMP recoil rate in an experimental detector. We will describe the 
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Standard Halo Model (SHM) and examine the behavior of the modulation signals in this 
model in Section IIII1 We will then examine how a dark matter stream affects this signal, 
first with the Sagittarius (Sgr) stream (Section II V|) . then generalizing to other streams 
(Section EJ. In Section |Vl we also briefly outline how various stream parameters can be 
extracted from a modulation signal. Our results are summarized in Section IVll 



II. DARK MATTER DETECTION 

WIMP direct detection experiments seek to measure the energy deposited when a WIMP 
interacts with a nucleus in the detector [41]. If a WIMP of mass m scatters elastically 
from a nucleus of mass M, it will deposit a recoil energy E = {p 2 v 2 /M){1 — cos#), where 
H = mMj (m + M) is the reduced mass, v is the speed of the WIMP relative to the nucleus, 
and 9 is the scattering angle in the center of mass frame. The differential recoil rate per 
unit detector mass for a WIMP mass m, typically given in units of counts/kg/day/keV, can 
be written as: 

„, d-R a(q) . . 



u>v min u 



where q = y2ME is the nucleus recoil momentum, a(q) is the WIMP-nucleus cross-section, 
p is the local WIMP density, and information about the WIMP velocity distribution is 
encoded into the mean inverse speed rj{E,t), 

r){E,t)= [ ^^-d 3 u. (2) 

Here 

% 

represents the minimum WIMP velocity that can result in a recoil energy E and f(u,t) 
is the (time-dependent) distribution of WIMP velocities u relative to the detector. For 
distinguishable WIMP populations (such as multiple streams), pr\ — > J2iPiVi i* 1 Eqn. (Q), 
where pi and rji are the local density and mean inverse speed, respectively, of some WIMP 
population indexed by i. 

In this paper, the use of the term "rate" will refer to the differential rate 1Z. Dates will 
be given as fractions of a calendar year {i.e. from January 1). For example, June 1 will be 
given as t = 0.415 years. For illustrative purposes, we will take a p = 1CT 42 cm 2 (defined 
in the following section) and use germanium as our detector element. In addition, as an 
example we take the WIMP mass to be m = 60 GeV and will use a local density of p = 0.3 
GeV/cm 3 for the smooth WIMP halo (streams will be in addition to the smooth halo, so 
the total local density can exceed this value). 



A. Cross-Section 

The cross-section for WIMP interactions is given by 



a{q) = a E 2 {q), 



(4) 
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where Cq is the zero-momentum WIMP- nuclear cross-section and F(q) is the nuclear form 
factor, normalized to -F(O) = 1; a description of these form factors may be found in Refs. 1*~12L 
4^ 1 . We will assume a purely scalar interaction, although the true interaction is likely to have 



both scalar and spin- dependent components. The results of this paper are not qualitatively 
affected by the choice of interaction: the primary difference would be a change in the recoil 
rates by an overall factor; the shape of the modulation would not change. 
For purely scalar interactions, 

4/y 2 

a = -^[Zf p + (A-Z)f n } 2 . (5) 

7T 

Here Z is the number of protons, A — Z is the number of neutrons, and f p and f n are the 
WIMP couplings to the proton and nucleon, respectively. In most instances, f n ~ f p ; the 
WIMP-nucleus cross-section can then be given in terms of the WIMP-proton cross-section 
as a result of Eqn. (jSJ): 

M A2 



where the \i v is the WIMP-proton reduced mass, and A is the atomic mass of the target 
nucleus. Again, for illustrative purposes, we have chosen a v = 10~ 42 cm 2 , a germanium 
(Z = 32, A w 73) detector, and a WIMP mass of m = 60 GeV. 



B. Velocity Distribution 

We will be examining two sets of WIMP populations: (1) a smooth galactic halo compo- 
nent, see Section Hill and (2) streams. We study the Sagittarius stream for illustration in 
Section HVl and generalize to all streams in Section[V] For both halo and stream components, 
we will use a Maxwellian distribution, characterized by a velocity dispersion a v , to describe 
the WIMP speeds, and we will allow for the distribution to be truncated at some escape 
velocity v esc , 

{_J_ ( _3A 3/2 p -3v 2 /2a2 r I , 
n csc \2ko$) 6 , ror |v| < v esc ^ 

0, otherwise. 

Here 

AU = erf(^)-2zexp(-z 2 )/7r 1 / 2 , (8) 
with z = v csc /vq, is a normalization factor. The most probable speed, 

v = <j2/3a v , (9) 

is used to generate unitless parameters such as z. For distributions without an escape 
velocity (v csc -> oo), N esc = 1. 

The WIMP component (halo or stream) often exhibits a bulk motion relative to us, so 
that 

/(u) = /(v obs + u) , (10) 

where v Q b s is the motion of the observer relative to the rest frame of the WIMP component 
described by Eqn. (JJJ); this motion will be discussed in the following section. For such a 
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velocity distribution, the mean inverse speed, Eqn. (J2J), becomes 



v (E,t) 



where 



i 



2N csc voy 
1 

2N csc v y 

o, 



eif(x+y) — erf(x — y) 



ye 



erf(z) — erf(x — y) — ^(y+z — x)e 



for z < y, x < \y—z\ 
for z > y, x < \y—z\ 
for \y — z\ < x < y+z 
for y+z < x 



x = v min /v , (12) 

y = Vobs/vo, and z = v esc /vo. Here, we use the common notational convention of representing 
3- vectors in bold and the magnitude of a vector in the non-bold equivalent, e.g. t> b s = |v bs|- 

Isothermal (Standard) Halo Model. For WIMPs in the Milky Way halo, the most fre- 
quently employed WIMP velocity distribution is that of a simple non-rotating isothermal 
sphere ^6|, also referred to as the Standard Halo Model (SHM). Typical parameters of the 
Maxwellian distribution for our location in the Milky Way are oshm — 270 km/s and t> esc 
= 650 km/s, the latter being the speed necessary to escape the Milky Way (WIMPs with 
speeds in excess of this would have escaped the galaxy, hence the truncation of the distri- 
bution in Eqn. ((7|)). More complex models, allowin g fo r e.g. anisotropy and triaxiality, may 
better match the actual dark matter halo; see Ref. |4^] and references therein. Such models 
do not qualitatively effect the results of this paper, as any smooth halo will give the same 
general behavior as exhibited by the SHM (non-smooth components, e.g. clumps, will result 
in streaming WIMPs, the focus of this paper). Then, for simplicity, we will take the SHM 
to describe the halo. The SHM will be discussed in Section HTT1 

Sagittarius stream. The Sagittarius (Sgr) dwarf galaxy is being absorbed bv the Milky 
Way and has a tidally stripped tail passing through the disk very near to us EMU!. This 
tail might provide a stream of WIMPs observable in dark matter detectors. While we are 
interested in streams in general, we will use this Sgr stream to illustrate how various stream 
properties affect the the signal in dark matter detection. In most cases, we will assume a 
dispersion of <7s g r = 25 km/s and an untruncated Maxwellian (t> esc — > oo) for this stream. 
The Sgr stream will be examined in Section IIVI 

General streams. In addition to the Sgr stream, there may be tidal streams from other 
dwarf galaxies being absorbed by the Milky Way or streams arising from late infall of dark 
matter. The dispersion cr str can vary, but will generally be much smaller than the observer's 
velocity t> b s in Eqn. (|10j). The Maxwellian for these streams will also be untruncated 
(f esc — * oo). General streams will be explored in Section fVl 



C. Motion of the Earth 

Due to the motion of the Earth around the Sun, v Q b s is time dependent: v Q b s = v b s (t). 
We write this in terms of the Earth's velocity V® relative to the Sun as 

Vobs(^) = v + V e [§i cosu(t - ti) + e 2 sin^(t - ti)] , (13) 

where u = 27r/year, v is the Sun's motion relative to the WIMP component's rest frame, 
V© = 29.8 km/s is the Earth's orbital speed, and E\ and £2 are the directions of the Earth's 
velocity at times t\ and ti + 0.25 years, respectively. Henceforth, as mentioned previously, 
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all times will be given as fractions of a calendar year (i.e. from January 1). With this 
form, we have neglected the ellipticity of the Earth's orbit, although the ellipticity is small 
and, if accounted for, would give only negligible changes in the results of this paper (see 
Refs. H EJ for more detailed expressions). For clarity, we have used explicit velocity 



vectors rather than the position vectors ei and e 2 used in Refs. [36, 40] and elsewhere (the 
position vectors are more easily generalized to an elliptical orbit); the two bases are related 
by k\ = — e 2 and — ei. 

In Galactic coordinates, where x is the direction to the Galactic Center, y the direction 
of disk rotation, and z the North Galactic Pole, 

e x = (0.9931,0.1170,-0.01032), (14) 
e 2 = (-0.0670,0.4927,-0.8676), (15) 

where we have taken e\ and e 2 to be the direction of the Earth's motion at the Spring 
equinox and Summer solstice, respectively, with t\ = 0.218 the fraction of the year before 
the Spring equinox (March 21). 

The time dependence of t> b s simplifies to the form 

v ohs (t) = ^Jv^+V£ + 2bv & V (B cosuj(t - t c ) , (16) 



where b = \fb\-\-b\ for b{ = e, • v , and t c is the solution of 

cosu(t c — ti) = — , smuj(t c — ti) = — . (17) 

VW+M VbJ+M 

The parameter b is a geometrical factor relating to the orthogonality of v with the Earth's 
orbital plane: b = sinA Q , where A Q is the angle between v Q and the normal to the orbital 
plane, with a maximum value of 1 when v is in the orbital plane and a minimum value of 
when v Q is completely orthogonal to the plane. We define the characteristic time as 

t c = the time of year at which t> b s is maximized. 

In the typical case in which the Earth's orbital speed is significantly smaller than the net 
motion of the WIMP population (or < u ) so that relative changes in f bsW are small, 
as with the SHM and (observable) streams, 



^obs(^) « f 



1 + b— costo(t — t c 

Vq 



(18) 



Standard Halo Model. Unlike the Galactic disk (along with the Sun), the halo has essen- 
tially no rotation; the motion of the Sun relative to this stationary halo is 

V©,SHM = V L SR + V 0jpec , (19) 

where vlsr = (0,220,0) km/s is the motion of the Local Standard of Rest and v Qi p ec = 
(10,13,7) km/s is the Sun's peculiar velocity. The characteristic time is t Cj sHM = 0.415 
(June 1) and the geometrical parameter b has the value of 0.49. The SHM will be discussed 
in Section Hill 
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Sagittarius stream. The Sgr stream is moving at approximately 300 ± 90 km/s relative 
to the galactic rest frame with a 3-velocity in this frame of 

v Sgr = 300 km/s x (0, 0.233, -0.970) . (20) 

The derivation of this velocity, described in Refs. 0,H(j, allows for significant uncertainties; 
however, as we will primarily use this stream simply as an example of streams in general, 
we will ignore these uncertainties and use only the stated values. The velocity of the stream 
relative to the Sun is then 

v 0,Sgr = V L SR + V 0ipec — V Sgr , (21) 

with vlsr and v 0jPec as given above, yielding t> 0i s g r = 340 km/s. The characteristic time for 
this WIMP population is then t CiSgr = 0.991 (Dec 28); b in Eqn. (O is equal to 0.53. The 
Sgr stream will be examined in Section IIVI 

General streams. For a general stream, we allow both the direction and magnitude of the 
velocity v str to vary. The characteristic time i C)Str can take any date and b can take any 
value from to 1; both parameters are completely specified by the stream direction v str . 
General streams will be explored in Section IS 



D. Annual Modulation 

It is well known that the count rate in WIMP detectors will exp erience an annual mod- 
ulation as a result of the motion of the Earth around the Sun |l5l Il6l|. In some cases, but 
not all, the count rate (Eqn. ((TJ)) has an approximate time dependence 



K(t) « K(t c + 0.25) 



1 +0(l)b— cosu(t 



(22) 



where the 0(1) factor may be positive or negative. For a cosine-like modulation, t c + 0.25 
years corresponds to the time of year at which the rate is approximately average; again, times 
are given as fractions of a year. Since we shall show that, in many cases, the modulation is 
decidedly not cosine-like, we divide the recoil rate TZ(t) more generally into a time-averaged 
component TZq and time-residual (modulation) component lZ m (t): 

K(t)=K +'R m (t)=Tlo[l + Xm(t)], (23) 

where Xmif) = U ^f > is the fractional (relative) modulation amplitude; by definition, 1Z m {t) 
averages to zero over one year. 

Typical WIMP recoil energies (on the order of 10's of keV) are comparable to that due 
to various background sources; experimental detectors thus have the difficult problem of 
needing to distinguish between background events and actual WIMP recoils, where the 
number of background events is much larger than the expected number of WIMP events. 
Presumably, however, such background does not exhibit the same time dependent behavior 
that is expected of WIMPs. For large detectors that can obtain high statistics, the WIMP 
modulation may be detectable even without precisely knowing the backgrounds present in 
the data; in this case, lZ m is therefore determined without knowing TZq. DAMA used a 
large Nal detector to search for such a signal and found an annual modulation This 
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signal is incompatible with the null results of other experiments, such as CDMS, under 
the conventional assumptions of elastic, spin-independent interactions and an isothermal 
halo 0] (with some exceptions These results, however, may be compatible for spin- 

dependent |47l liij l or inelastic interactions (jjl. 5oL l5ll ] . Interpretation of the DAMA signal 
for non-standard halos is discussed in Refs. [5_2|, |5_3j ; DAMA also examined the possibility of 
streams in Ref. The presence of an energy cutoff due to a stream in this modulation 
signal would provide a "smoking gun" for WIMPs. 

The primary purpose of this paper is to examine the annual modulation due to WIMP 
streams and constrast this modulation with that due to the smooth galactic component. 
We will begin by examining the smooth galactic halo component in Section IHII We will 
then compare how a stream modulation differs from the galactic component and study 
how various stream parameters affect the modulation, first using the Sagittarius stream for 
illustrative purposes in Section [JY] and then with more general streams in Section ["""J 



III. ISOTHERMAL (STANDARD) HALO MODEL 

As noted in the previous section, the most frequently employed WIMP velocity distribu- 
tion for the Milky Way halo is that of a simple isothermal sphere also referred to as 
the Standard Halo Model (SHM), given by Eqn. (J7j) with oshm = 270 km/s and v esc = 650 
km/s. For these parameters, previously discussed in Section Hi Bl and the parameters given 
in Section 111 CI (we will also take the local WIMP density for this smooth halo component 
alone to be pshm = 0.3 GeV/cm 3 in this and all following sections), the modulation of the 
rate is well approximated by a cosine over a large range of energies, with a relative residual 
rate (see Eqn. (HHJ)) 

f-0.034(l-4), forx<a;„, 
Xm(x,t) « cosco(t-t c ) x <^ x p (24) 

1 ' \0.014(J-l)(^ + 3.7), forx p <x<z, 1 1 

where t c = 0.415 (June 1), z = 2.95, and x v = 0.89 is the value of x at which the phase of 
the modulation reverses (recall x = v m - m /v with t> min oc VE). For x>z, the relative size 
of the modulation grows and the shape becomes non- cosine-like; however, the overall rate 
becomes highly suppressed. From Eqn. (J2H), we can see two features: 

1. The modulation is only a few percent of the average count rate. Thus, a large number 
of events are required to observe a modulation of the rate in a detector. 

2. For x < x p (low recoil energies), the cosine is multiplied by a negative factor: the rate 
is minimized at a time t c . For x > x p (high recoil energies), the reverse is true: the 
rate is maximized at time t c . See Ref. [55] for a discussion of this phase reversal. 

The detector type (composed of an element with mass M) and WIMP mass m arise in 
Eqn. ([2*4]) only through a constant scaling of the parameter x = with the square root of 
the recoil energy E via Eqn. (J"JJ). Thus, the modulation behavior of Eqn. ([2*4]) and the two 
features listed above are expected to be observed in any detector, for any WIMP mass, but 
at different energy scales. For the sake of illustration, we will assume a WIMP mass m = 60 
GeV and a Germanium (A m 73) detector. In this case x = 0.25^ E/lkeV; the parameter 
value of x p = 0.89 at which the phase of the modulation reverses in Eqn. ()24j) corresponds to 
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FIG. 1: The modulation due to the Standard Halo Model (SHM) and Sagittarius (Sgr) stream 
separately, as well as the combined modulation, for a Germanium detector. Each component is 
shown at a recoil energies of 10 keV and 39 keV. The WIMP mass is assumed to be 60 GeV and 
the Sgr density is taken to be 5% of the SHM density. 



an energy E, p = 13 keV. Then for E < 13 keV, there will be a cosine-like modulation in the 
recoil rate at the few percent level, minimized around June 1; see Figure ^ for an example. 
For 13 keV < E < 110 keV, a cosine-like modulation will occur, but will be maximized 
around June 1; again see Figure H For 22 > 110 keV, the overall rate becomes insignificantly 
small (compared to lower energies). Experiments such as CDMS 0| and EDELWEISS p, 
which use Germanium detectors, have thresholds low enough to potentially see this phase 
reversal, but are not expected to observe enough events to discern the small (few percent] 
modulation effect. Experiments using other elements in their detectors, such as ZEPLIN jy] 
and XENON 0] (using Xenon), likewise could see such a reversal if the exposure is sufficient 
to detect the modulation effect. 

Both features, the small modulation amplitude and phase change, arise predominantly 
from the mean inverse speed factor 77 in Eqn. (JTJ; rj is shown in Figure |2] at both t c and t c 
+ 0.5 years. At t c , the Earth is moving the fastest relative to the SHM, shifting the WIMP 
velocity distribution, Eqn. (jTDjl. to higher velocities compared to t c + 0.5 years (when the 
Earth is moving slowest). The total recoil rate (over all energies) is highest at t c , but 
the shift to higher velocities depletes the number of low velocity WIMPs and low energy 
recoils. However, the small change in our velocity throughout the year relative to our velocity 
through the SHM, along with the large velocity dispersion of WIMPs (comparable to our 
motion through the halo), yield only small changes in 77; the fractional change is of order 
^0~shm ^ or a ^ ener gi es - The large velocity dispersion reduces the modulation effect because, 
even in our frame, there are significant numbers of WIMPs incident on Earth from all 
directions (but somewhat more likely from the "forward" direction): in some directions, the 
WIMP velocity is actually reduced at t c and increased at t c + 0.5 years. The geometrical 
factor b in Eqn. ()16|) (relating to the angle of the Sun's motion through the galaxy with 
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FIG. 2: The mean inverse speed rj, given by Eqn. (J2J), for the SHM and Sgr stream at t c and 
t c + 0.5 years, where t c is different for the two WIMP populations. The characteristic time t c is 
the time of year at which the Earth is moving fastest relative to the given WIMP population (see 
Section Til C|) . For the SHM, rj at t c ,SHM (June 1) is slightly larger at higher recoil energies and 
smaller at lower recoil energies than at £ C ,SHM + 0.5 years (Dec 1); the change in the phase occurs 
at an energy -E^shm = 13 keV. In the SHM, the fractional change in r\ is of order „ ^ M ( onr y a 
few percent) for all energies. For the Sgr stream at t Ci s gr (Dec 28), rj is flat up to a cutoff near 
the characteristic energy E c a, gI = 39 keV (see Eqn. (|28j0 ; the cutoff is softened by the velocity 
dispersion of WIMPs in the stream. At i Ci s g r + 0.5 years (June 29), the cutoff energy has decreased 
to around 35 keV due to the lower maximum WIMP velocities, but the plateau height at energies 
below this has increased relative to t c ,sgr- For rj at energies below the cutoff, the relative change 
over the year is also of order — - — . However, due to the shift in the cutoff, the change in n is 



not suppressed by this velocity factor over a small range of recoil energies near 39 keV; the change 
here is 0(1) relative to the average at that energy. 

the Earth's orbital plane), here having the value 0.49, further suppresses the modulation 
amplitude. 

As is apparent from Eqn. (|24|) . the modulation is symmetric in time about t c (June 1) 
for the SHM. The rate is always at an extremum on this date, with the other extremum 0.5 
years later (Dec 1), as seen in Figure El We shall see in the next section that there is no 
symmetry when additional WIMP populations, such as streams, are present in the halo: the 
date of the rate extrema changes with energy and the minimum and maximum rate need 
not occur 0.5 years apart. 



"0.SHM 



IV. 



SAGITTARIUS STREAM 



Beyond the smooth SHM dark matter distribution, the galaxy contains additional struc- 
ture, such as streams arising from late infall of dark matter or tidal disruption of dwarf 
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FIG. 3: The dates at which the minimum and maximum recoil rates occur for various recoil 
energies. For the SHM alone, as demonstrated in Figure Q the modulation is symmetric about 
£c,SHM (June 1), with a minimum at this time for energies below E Pi shm. = 13 keV and a maximum 
here for energies above this; the other extremum is always at i c ,SHM + 0.5 years (Dec 1). The Sgr 
stream component alone is symmetric about £ c ,Sgr (Dec 28); but when the SHM and Sgr stream are 
combined (shown here with a 5% stream density relative to the SHM density), the total modulation 
becomes asymmetric. The asymmetry is apparent at 30-50 keV, where the maxima and minima 
are not 0.5 years apart. 



galaxies being absorbed by the Milky Way. While we are interested in examining streams 
in general, in this section we will use the Sagittarius (Sgr) stream 2^, 23|, 3 as an example 



to illustrate how various stream parameters affect the modulation signal in a detector. As 
discussed previously, the Sgr stream is a stream of WIMPs associated with a tidally stripped 
tail of the Sgr dwarf galaxy that is currently being absorbed by our galaxy. This tidal tail 
passes through the Milky Way's disk very near to us, so the WIMPs in the Sgr stream are 
potentially observable in direct detection experiments. Detection of WIMPs from the Sgr 
stream has been discussed in Refs. |39j, |40j ; here, we expand upon those discussions. We will 
discuss the dependence of the annual modulation in the count rate on the recoil energy, the 
binning of the recoil energy, the stream density, and the stream dispersion. In Section we 
will examine how the results illustrated with the Sgr stream can change with more general 
streams. From the parameters discussed in Section III Bl and Section III CI the characteristic 
time for this WIMP population, i.e. the date at which we are moving fastest relative to the 
stream, is t c = 0.991 (Dec 28). The geometric parameter b in Eqn. (|16|) . associated with the 
angle of the stream relative to the Earth's orbital plane, is equal to 0.53. 

We first examine the basic behavior of a stream by neglecting the velocity dispersion 
(cr„ = 0). In this case, the mean inverse speed r], Eqn. (j2J), has the constant value - up 
to a cutoff energy E co (t), 

r] (E,t) = — i --r6[E co (t)-E) , (25) 

Vobs (Z J 
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where 6 is the Heaviside step function. The cutoff energy corresponds to t> min = v ^ s (t) ( see 
Eqn. ©), so 

E co (t) = E C [1 + A c cos u(t - t c )\ , (26) 

where 



2&^0 ,SgrVe 

; 0,Sgr 

We define a characteristic energy 



A c = 7 "Tt/2 - ( 27 ) 



£ c - ^(Kbs(t)] 2 ) = ^ « Sgr + Vi) . (28) 

The relevant velocities have been defined in Section 111 CI fsee. e.g., Eqns. (j2Ti|) and (j2U))- 
The characteristic energy E c is also the average cutoff energy for the case of zero velocity 
dispersion discussed here, but the definition in Eqn. (J2~%j) may still be used when a velocity 
dispersion is included (a v ^ 0; see below), in which case there is no hard cutoff energy at 
E co (t). While this definition can also be used to define an E c for any WIMP population, E c 
is not a useful quantity for a WIMP population with a large velocity dispersion (a v >v Q ), 
such as the SHM, as there is no associated rapid drop off in the count rates near that energy 
such as with a stream (which has a v <C v Q ). Hence, we will only use the characteristic energy 
of streams in our discussion. We wish to emphasize the importance of the step function in 
Eqn. (J25j) : the presence of the step and the fact that its position varies with time have key 
consequences that will be seen in the discussion that follows. 

As with the SHM modulation discussed in the previous section, the modulation from the 
stream (arising from Eqn. ()25j0 occurs for any detector type, composed of an element with 
mass M, and any WIMP mass m. However, the characteristic energy E c given by Eqn. ()28|) 
and, thus, the location of the step as given by Eqn. (j2T)j) depend explicitly on M and m 
(recall /i = Mm/(M + m) is the reduced mass). The choice of M and m do not change 
the qualitative behavior of the modulation, only the energy scales at which various effects 
occur. For illustrative purposes, we take a 60 GeV WIMP and a Germanium detector, for 
which E c = 39 keV for the Sgr stream; our results will apply to other detectors and WIMP 
masses, just at different energies. 

For the dispersionless case, 

Xm (E,t) « -0.046 cos uj(t - t c ) , for E < E coMn (29) 

and there are no recoils for E > E co ^ maiX , where -E C o,min and -E C o,max are the minimum and 
maximum cutoff energies, occurring at t c + 0.5 years and t c , respectively. Below -E co , m in, the 
relative variation in the recoil rate is of order — — and the modulation of the rate is cosine- 
like with a minimum at t c ; this behavior is similar to that seen in the SHM for energies 
below the phase reversal energy E p . For E > E co (t), the WIMPs in the stream are not 
moving sufficiently fast to produce a recoil of energy E. However, for -E co , m in < E < E COjmax , 
there are times during the year at which E > E co (t) (71 = 0) and there are times at which 
E < E co (t) (71 > 0). In this case, the size of the modulation is 0(1) relative to the average 
recoil rate due to the step in Eqn. (l2*oT) . much larger than the 0( v ® ) effect arising from 
the time dependence in — Krr that is apparent in Eqn. (1181) : the modulation in this energy 

^obs \t) 

range can be quite large and non- cosine-like. 

The behavior of Eqn. (|2*3j) approximately holds for non- zero velocity dispersion a v as well, 
provided the dispersion is not significantly larger than the changes in f b s (£) («'-e. a v <bV e ), 
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although the cutoff softens for non-zero a v . For illustrative purposes, we take the velocity 
dispersion of Eqn. (JJJ) to be a v = 25 km/s for the Sgr stream (corresponding to Vq = 20 
km/s from Eqn. Q), although we will examine variations in this parameter in Section Ft V El 
For this case, r] is shown in Figure|2]at both t c (Dec 28) and t c + 0.5 years (June 29). In the 
figure, the softened step function is apparent, with a cutoff around E c = 39 keV. The height 
of the step and the cutoff energy vary with time. At t c , the Earth is moving the fastest 
relative to the stream. This leads to a larger range of recoil energies as opposed to t c + 0.5 
years, when the Earth is moving the slowest, so the step occurs at a larger energy at t c (~43 
keV) than half a year later (~35 keV). The higher WIMP velocities, by spreading the recoils 
over a larger energy range, lead to less recoils at any given energy and a lower step height. 
The relative variation in the step height is of order b v<s (~5%), similar to the relative 

u Q,Sgr v n 

variation in rj for the SHM (also a few percent). However, the variation around the step is 
unlike any seen in the SHM: rj goes from nearly zero at t c + 0.5 years to essentially the full 
step height at t c . Just as with the dispersionless case previously discussed, the dominant 
contribution to the modulation in this energy range is the shift on and off the step, not 
the variation in the step height itself. Thus, the relative variation in the recoil rate (recall 
1Z{t) oc rj(t)) is not suppressed by the above velocity factor and Xm ~ C(r)- This is a large 
amplification in the signal and could potentially yield an observable signal for the Sgr stream 
at a comparable level to the SHM, even if the Sgr stream is significantly less dense. 

To examine the recoil rate for the Sgr stream and compare to the SHM, we must assume 
a local density psgr for the stream. For illustrative purposes, we will assume a stream 
density 5% that of the SHM (p Sgr = 0.05pshm)- We note this is optimistic and that the 
local Sgr density is likely at most a few percent j4£j. However, 5% is not unreasonable for 
other structure in the dark matter distribution, such as clumps (3?J. The various effects we 
examine still mainly apply for lower densities, just to smaller degrees. We will, however, 
examine variations in the density in Section TlV Dl 

We will take psgr = 0.05psHM and a v = 25 km/s to be our fiducial values in Sections llV A[ 
IIVCI and then we will vary them in Sections IIVDI & IIV El respectively. We will examine 
the dependence of the modulation in the count rate as a function of the recoil energy, the 
binning of the recoil energy, the stream density, and the stream dispersion. We wish to 
reiterate the following definitions, which play an essential role in the remaining discussions: 

the characteristic time t c is the time of year at which Earth is moving fastest rel- 
ative to some WIMP population, formally defined by Eqn. (fTTj). and 

the characteristic energy E c is the approximate energy at which a rapid drop in 
the recoil rate occurs for streams, formally defined by Eqn. (j2EJ). 

We also reiterate that, for illustrative purposes, we take a 60 GeV WIMP and a Germa- 
nium detector, for which E c = 39 keV for the Sgr stream; our results will apply to other 
detectors and WIMP masses, just at different energies. 

A. Contrasting and Summing the Halo and Stream Modulations 

Here we will compare the modulation due to a stream, using the Sgr stream as an example, 
with the modulation due to the SHM as discussed in Section IIHI and we will examine the 
total modulation (sum of the SHM and stream components) that a detector will observe. 
Figures demonstrate out results. 
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For a 5% Sgr stream with a v = 25 km/s, the residual rate (modulation) is shown in 
Figure ^ at recoil energies of 10 keV and 39 keV. At 10 keV, the modulation is essentially 
a cosine minimized on December 28 (t c ,sgr) an d maximized on June 29. This behavior is 
similar to the low energy behavior of the SHM, although the SHM has a different t c (June 1) 
and, hence, does not peak at the same time. At 39 keV, near the characteristic energy E C} s gT 
for the Sgr stream, we see the modulation becomes extremely large for the stream, with an 
amplitude nearly four times larger than at 10 keV, even though the nuclear form factor in 
Eqn. (jlj) generally causes the recoil rate to decrease at higher energies. This amplification in 
the amplitude is due to the r\ step edge crossing discussed previously. The phase at 39 keV 
is reversed from that at 10 keV: the maximum now occurs on December 28 (£ c ,sgr) rather 
than on June 29. However, unlike the SHM, the stream modulation is not cosine-like at all 
energies. While this particular instance may look somewhat cosine-like, most energies in 
the cutoff range around E c are decidedly not cosine-like; a smaller o v would also make the 
modulation at E c $ gx l ess cosinusoidal (we shall see this in Section II V Kj) . 

In general, even if the WIMP population components (e.g. SHM and Sgr stream) indi- 
vidually produce cosinusoidal modulations, the total modulation seen in detectors is not 
cosinusoidal. At 10 keV, the total modulation shown in Figure ^ does look somewhat cos- 
inusoidal; this, however, is a coincidence due to the nearly six month difference between 
the t c 's of the two components. At E C) s gT = 39 keV, the total modulation is clearly not 
time-symmetric. 

While the SHM and Sgr stream always have their maxima or minima at their respective 
t c 's for any energy (but possibly changing between the maximum and minimum at some 
specific energy), the combined modulation has extrema occurring on dates that vary with 
energy. In Figure El the SHM can be seen to always have the extrema on June 1 (£ c ,shm) 
and December 1, but with the extrema reversing at an energy E P) shm = 13 keV. With a 5% 
stream in addition to the SHM, the combined modulation is mostly the same as for the SHM 
alone except near -E Pi shm and -E c ,Sgr- At low energies (E<8 keV), the SHM has the most 
significant contribution to the modulation and the extrema occur on nearly the same dates 
as the SHM alone. At intermediate energies (15 keV <E< 30 keV), the SHM is again the 
most significant component and the extrema of the total modulation are near that of the 
SHM's. Above 45 keV (above the maximum cutoff energy -E co ,max for the Sgr stream), there 
is essentially no contribution from the stream, so the extrema match the SHM. In the regions 
near E Pi sbm and E Cj s gr , however, the dates at which the maximum and minimum count rates 
occur are different for the combined SHM & Sgr stream modulation than they would be for 
the SHM alone. The SHM modulation amplitude varies smoothly across zero at -E Pi shm and 
is small at nearby energies, as can be seen in Figure El The Sgr stream modulation amplitude 
does not vary significantly at these energies, so the stream contributes significantly to the 
modulation and leads to a smooth variation in the extrema dates. Around and between 
E co ,mm and E COimax for the stream (33-45 keV), the total modulation fluctuates wildly due to 
the cutoff in rj. The modulation here is noticeably asymmetric in Figure El as the extrema 
are not 0.5 years apart. 

As noted above, several properties can emerge from multiple components that are not 
present in the individual components: 

• The combined modulation is generally non-cosinusoidal, even if the modulation of each 
individual component is; 

• The combined modulation need not be time-symmetric, even if each individual com- 
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FIG. 4: Dependence on recoil energy. The solid (blue) line shows the modulation due to only the 
SHM at a recoil energy of 39 keV (f); the remaining lines include a 5% Sgr stream in addition to 
the SHM at various recoil energies. The shape of the modulation is dependent on the recoil energy 
and clearly exhibits non- cosine-like behavior over a range of energies. The contribution from the 
stream falls off at energies above 39 keV, the characteristic energy -E c ,Sgr for the Sgr stream (see 
Figure EI). 



ponents is; and 

• The minimum and maximum recoil rates need not occur 0.5 years apart. 

These three effects, if observed, are potentially evidence for some structure (streams) in 
the local WIMP population. These effects will be apparent in the following sections, as we 
examine the behavior of the modulation for different recoil energies, recoil energy binning, 
stream density, and stream dispersion. For these sections, we will be examining only the 
combined modulation of the SHM and Sgr stream (not individually). 



B. Recoil Energy 

Here we discuss the dependence of the total modulation on the recoil energy. As seen in 
Figure HJ the shape of the modulation is highly dependent on this energy. In this figure, we 
have chosen to show the total modulation at representative energies below, at, and above 
39 keV, the characteristic energy of the Sgr stream for 60 GeV WIMPs in a Germanium 
detector. For comparison, the cosine-like contribution from the SHM alone is shown at 39 
keV (solid/blue line). 

Well below the characteristic energy of the stream, both the SHM and stream modulations 
are cosine-like, although peaking at different times. In FigureEJ we see that the combination 
is approximately cosine-like at 29 keV, a consequence of the two separate modulations being 
nearly in phase with each other (both the SHM and Sgr stream peak in the summer at this 
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energy). For the SHM (and other WIMP populations), the cosine-like modulation comes 
from the expansion of the SHM's speed relative to the Earth Vob s (t) in the powers of — <C 1 

y'2 

as given in Eqn. (TTH1) . This equation has corrections of 0(-§-). Although not noticeable in 

^0 

the figure, the combined SHM + Sgr stream modulation at 29 keV differs from a true cosine 

v 2 

by a larger amount than these 0(-§-) corrections. A stream with a different t c would yield 

a modulation that does not appear as cosine-like. 

At 39 keV (E Cy s gr ), the 0(1) relative modulation of the stream leads to a contribution 
to the modulation by the stream comparable in magnitude to the contribution from the 
SHM, shifting the peak date by nearly three months and greatly decreasing the modulation 
amplitude relative to the SHM dominated behavior at 29 keV or SHM component of the 
modulation alone at 39 keV. This is a highly significant result: even though the stream is 
much less dense than the SHM, its modulation is just as large as that of the SHM here. The 
stream's effect on the total modulation is likely to be seen here well before other energies. 
The asymmetry in the modulation is apparent at this energy. The contribution to the total 
modulation from the stream is also significant at energies near to E C} s gr , as can be seen at 
34 keV and 44 keV (the latter being more apparent near the minima in the figure). 

At 49 keV, above -Eco.max for the stream, there is essentially no contribution from the 
stream and the total modulation assumes the cosine-like form of the SHM. The amplitude 
is reduced from the case at 29 keV, which had contributions from the stream, and is also 
slightly smaller than the SHM component alone at 39 keV, due to the form factor in Eqn. (0J) 
that lowers the count rate at higher energies. 

A stream can have a significant effect on the modulation near the characteristic energy E c 
even if the stream is far less dense than the smooth background halo distribution (here, the 
SHM). The stream can significantly change the peak date and the shape of the modulation 
(becoming quite non-cosinusoidal) near this energy. Away from this energy, however, a low 
density stream has only a mild impact, if any, on the modulation. The question of how well 
experiments are able to distinguish the modulation behavior near E c , due to the limited 
energy resolution in a detector, is addressed in the following section. 



We examine how the modulation appears when averaged over some range of recoil ener- 
gies, as would be observed when experimental data is binned, and demonstrate the results 
in Figure We have seen that the Sgr stream has a significant effect on the modulation 
at least over the recoil energy range 34-44 keV (see Figure EJ) . As detectors have limited 
energy resolution, experimental results are often binned in recoil energies. If the bins are 
much larger than the limited width of the stream's 0(1) relative modulation effect, it would 
be difficult to study this behavior or even observe it. In Figure the modulation from the 
SHM and a 5% Sgr stream is shown for multiple bin sizes (AE), where the bins are centered 
at 39 keV, the characteristic energy E Ct s gT (corresponding the energy cutoff; see Figure EJ). 
The average recoil rate, 



C. Recoil Energy Range (Binning) 




(30) 



is used in each case. For comparison, the SHM component alone is shown at 39 keV. 
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FIG. 5: Dependence on binning of recoil energy. The average modulation observed over an energy 
range centered at 39 keV (E C} Sg r ) with a width AE. The solid (blue) line is for the SHM only with 
zero bin width (f); the remaining lines include a 5% Sgr stream in addition to the SHM. Only 
for large bins, AE > 10 keV, does the shape of the modulation differ significantly from that of the 
differential recoil rate (equivalent to zero bin width: AE = keV). 



An infinite resolution detector, with AE = 0, will observe the actual recoil rate at E c ^ gr - 
For bin widths of 2 & 5 keV, there is very little deterioration in the modulation signal. Even 
for a 10 keV bin, the modulation signal is similar to the AE = case. For AE = 20 keV, 
the signal is highly deteriorated, but the modulation observed is still very different from the 
SHM case. The current generation of Ge detectors have energy resolutions on the order of 
1-2 keV, so it is clear detector resolution will not inhibit observation of the stream's effect 
on the modulation around E C) s gT . 



D. Stream Density 

The extent to which a stream affects the annual modulation depends upon the density of 
the stream. The modulation due to the SHM + Sgr stream is shown in Figure EH for several 
different stream densities, all at the characteristic energy -E Cj s g r- For comparison, the SHM 
only case is shown by the solid (blue) line. At this energy, the modulation of the SHM 
peaks on June 1; that of the Sgr stream is minimized on June 29. The stream obviously 
significantly contributes to the overall modulation at this energy for the 5% relative density 
we have been using. 

While a 5% density may be reasonable for some streams j3?J (the caustic ring model 
of P. Sikivie and collaborators predicts a stream as large as 75% of the local density; see 
Ref. [38| and references therein), it is optimistically large for the Sgr stream. Instead, a 
density of a few percent or lower may be more reasonable. Would such low densities be 
observable? In the figure, streams at 5% and 10% densities differ significantly from the 
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FIG. 6: Dependence on stream density. The modulation for a variety of stream densities at a 
recoil energy of 39 keV (E c ,sgr)', stream densities are given relative to the (fixed) SHM density. 
The solid (blue) line represents the SHM alone. Even at only 5% of the density, the stream has 
a comparable modulation amplitude to that of the SHM alone near the cutoff energy, yielding a 
total modulation quite different than would be expected for the SHM alone. 

SHM component. In these cases, the modulation is non-cosinusoidal and asymmetric and 
has a different amplitude and peak date than the SHM component. For the 5% case, the 
maximum rate occurs 3 months prior to that of the SHM component. For a 2% stream, the 
modulation shape and peak date do not differ significantly from that of the SHM. However, 
even though the stream density is only 1/50 that of the SHM, the amplitude of the total 
modulation is still decreased by more than 30% from that of the SHM alone near E c ^ gr . 
The implication then is, due to the 0(1) relative modulation of the stream near -E c ,sgr, a 
relatively small stream can still have a significant effect on the experimental results. 



We examine the dependence of the modulation on the velocity dispersion a v of the stream, 
illustrated in Figure Up to this point we have taken the dispersion to be 25 km/s, and we 
now allow it to vary. The SHM and Sgr stream have very different modulation behaviors: the 
SHM modulation is always very nearly cosinusoidal, with a small amplitude and a reversal 
of the modulation phase at an energy -E Pi shm- The stream modulation, on the other hand, 
has three different behaviors: it is cosinusoidal with a small amplitude at energies below the 
characteristic energy E c ^ gT , it is large (0(1) relative to the average rate) and non-cosinusoidal 
around -E Cj s g r, and it is non-existent above this energy. The differences between these two 
components are mainly due to the large difference in their velocity dispersions, leading to the 
following two consequences: (1) For a dispersion of a WIMP population significantly smaller 
than the net motion u b s of that population relative to a detector, or a v C %, the recoil 



E. Stream Dispersion 
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FIG. 7: Dependence on stream dispersion. The modulation at (-£7 c ,Sgr) is shown for various stream 
velocity dispersions (a v ). The solid (blue) line is for the SHM only; the remaining lines include a 5% 
Sgr stream in addition to the SHM. For a v = km/s, the stream has a step function contribution 
on top of the cosine-like modulation of the SHM. Increasing velocity dispersion softens the step 
function like behavior. 



energy spectrum develops a relatively rapid drop near the characteristic energy E c (true for 
the Sgr stream but not for the SHM). (2) For a dispersion on the order of the variation in 
v bs{t) due to the Earth's orbital motion, or cr v <bV<£, the variation of the location of this 
dropoff leads to relatively large, non-cosinusoidal modulations in the count rates near E c 
(again, this applies to the Sgr stream but not the SHM). This second condition essentially 
requires that the shift in the location of the step (~ bV®) be comparable or larger than the 
width of the step dropoff itself (~ a v ); if this condition is not satisfied, the modulation near 
E c remains a relatively small effect (only a few percent of the total rate, as occurs at other 
energies). While the presence of a dropoff is signified by the condition of (1), the (usually) 
stronger condition of (2) is necessary to make the dropoff rapid enough to observe the large 
modulation effect. 

For the stream, T](E) is similar to a step function, with the position of the step, Eqn. (J2BJ), 
shifting in time; see Figured The shift in the step is proportional to 5 V®, the change in 
stream velocity v h s (t) (recall < b < 1 is a geometrical factor dependent upon the stream 
direction). The edge of the step, however, is softened, with the fall off occurring over an 
energy range proportional to o~ v . As long as the velocity dispersion is not significantly 
larger than the changes in v b s (t) (i.e. (j v <bV®), the edge is sharp enough so that, at some 
energies, a large portion of the edge crosses the given energy throughout the year. Note how 
in Figure El f] is nearly atop the step at 39 keV on December 28 (t c ,Sgr) for the Sgr stream, 
but is beyond the step on June 29. When a v ^> 6V®, the step is washed out, as can be 
seen by r\ for the SHM in the same figure. The small dispersion that is a characteristic of 
streams leads to the 0(1) relative modulation (arising from the crossing of the step) that is 
not present in the SHM with its large dispersion. 
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The SHM + 5% Sgr stream modulation is shown in Figure [7| for multiple stream dis- 
persions at a recoil energy of 39 keV (E Ct s gr ). For comparison, the SHM only case is again 
shown by the solid (blue) line. The step like behavior in 77 for the stream is apparent as the 
dispersion goes to zero. At a v = 0, where r\ is a true step function, the step edge, E co , can 
be seen to have cross 39 keV in late March and early October. From March to October, 
the stream cutoff energy is below the recoil energy being observed, so that the stream does 
not contribute to the count rate (note how the shape of the modulation on these dates is 
identical to the SHM only component). From October to April, the cutoff energy is above 
the recoil energy being observed, so the stream also contributes to the count rate. 

As <j v increases, the sharp changes in the modulation in March and October are softened. 
o v = 12.5 km/s still looks somewhat similar to the a v = case. The a v = 25, 37.5, & 50 
km/s cases all have a v greater than &V® = 16 km/s; the abrupt changes in the recoil rate are 
now gone. However, a v is not significantly larger than 6V® and an appreciable portion of the 
rj step still crosses this recoil energy: the lighter stream still has an effect on the modulation 
comparable to that of the SHM, as evident by the large difference between these cases and 
the SHM only case. 



As noted by Refs. [37|, |56j, |57[, the stream's velocity dispersion is more likely to be 
anisotropic, in which case the distribution of Eqn. (JJJ) is not valid. However, including 
anisotropic models of the stream's velocity distribution would not significantly affect the 
results discussed in this paper as the velocity dispersion of the stream will still be significantly 
smaller than the stream velocity relative to the observer. In such models, the dropoff in the 
recoil rate would still be present, as well as a modulation of both the recoil rate and location 
of the dropoff. Including such models would simply lead to modest changes in the shape of 
the dropoff in the mean inverse speed (see Figure [2J) and, hence, the recoil rate around the 
characteristic energy E c . 



V. GENERAL STREAMS 

Having examined various stream parameters using the Sgr stream in the previous section, 
we wish to expand our results to streams in general. The Sgr stream modulation features 
previously discussed likewise occur for other streams: (1) The phase of the stream mod- 
ulation is associated with a characteristic time t c (typically the time of year at which the 
recoil rate is maximized or minimized) that is independent of the SHM; (2) There exists a 
characteristic energy E c associated with a rapid dropoff to zero in the recoil rate; (3) At 
energies less than the characteristic energy E c , the modulation is nearly cosinusoidal, with 
a minimum at t c and a (small) amplitude of 0(^-) with respect to the average rate; and 
(4) The modulation near E c is not cosinusoidal, is of 0(1) relative to the average rate, and 
can be comparable in amplitude to that of the overall halo modulation. While any stream 
will have these features (we are assuming a small velocity dispersion, or a v <V (S , as would 
be expected for typical streams), the characteristic time t c and energy E c depend upon the 
stream. Given a stream velocity v str , these values may be determined from Eqn. (|17|) and 
Eqn. respectively (using v &tStT in place of t^sgr m the latter equation). We will examine 
how the stream speed and direction affect the characteristic time and characteristic energy, 
and consequently the modulation. 
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Stream Speed (km/s) 

FIG. 8: Stream speed. The characteristic energy E c , defined by Eqn. (|'28[). as a function of stream 
speed. E c is approximately the highest possible recoil energy from a WIMP in the stream. E c is 
independent of the stream direction and depends only on the speed of the stream relative to the 
Sun (f0,str; solid/blue line). The 340 ± 70 km/s speed of the Sgr stream is shown (red dot with 
wide pink line indicating the la region), corresponding to an E c range of 25-55 keV with central 
value 39 keV. 

A. Stream Speed 

The characteristic energy E c of a stream, Eqn. (|2*Hj) with ^o,Sgr replaced with a more 
general stream speed v QjStT , is dependent only upon the speed of the stream relative to the 
Sun and is independent of the stream direction. For a stream moving much faster than 
Earth's orbital velocity (v QtStr 3> V®), the characteristic energy is proportional to the square 
of the stream speed relative to the Sun (E c oc Vq str ), so faster streams rapidly lead to higher 
E c 's; this behavior is demonstrated in Figure El 

The E c determined from a modulation signal can be used to determine the speed of 
the stream to an improved degree over that afforded by limited alternative observations. 
The 300 ± 90 km/s Sgr stream speed estimate relative to the galactic rest frame, discussed 
in Refs. 0, H(j , is not based upon direct observation of Sgr stellar material in the local 
neighborhood, but on extrapolating stellar stream observations in other areas of the Milky 
Way, leading to the relatively large uncertainty. The 300 ± 90 km/s speed relative to the 
galactic rest frame corresponds to a speed relative to the Sun of v Qt sgr = 340 ±70 and, hence, 
to a 25-55 keV range of E c , as shown in Figure |H1 This energy range is much larger than 
the energy resolution of a WIMP detector; current detector resolutions of ~1 keV would 
yield a stream speed accurate to about ±10 km/s (in the Sgr stream case). To derive the 
stream speed from E c , using Eqn. (J2~%j) (where M is the nuclear target mass and fi is the 
WIMP-nucleus reduced mass), the WIMP mass must already be reasonably well known. 



22 



T 




: i i i i i i 

100 200 300 

Stream Angle (degrees) 



FIG. 9: Stream direction. The characteristic time t c as a function of stream direction; t c is 
independent of the stream speed. Each curve shows the dependence on one of the three angles, 
in the xy-, xz-, and yz-planes (e.g., the angle 4> in the xy-plane corresponds to the direction 
cos (j)x + sin <f>y in Galactic coordinates). 

1. Determining the WIMP Mass 

Alternatively, if the WIMP mass is unknown, but a dropoff in the recoil rate can be 
associated with a known stream (known via other observations), the energy at which that 
dropoff occurs, E c , can be used to determine the WIMP mass. For instance, if we observe a 
dropoff at E c = 30 ± 2 keV in a Germanium detector and associate it with the Sgr stream 
(with speed 340 ±70 km/s), we would derive a WIMP mass of 48 ±22 GeV. If further stellar 
observations and modelling of the Sgr stream improve the speed uncertainties to ±20 km/s, 
the mass determination would improve to 48 ± 7 GeV. Certainly, the presence of a dropoff 
in the recoil rate, which does not require a modulation effect to be observed, can provide a 
useful tool for determining the WIMP mass. However, care must be taken in associating such 
a dropoff with a specific stream: observation of a 5% density stream through such a dropoff 
alone cannot be assumed to be due to the Sgr stream rather than some as yet unknown 
other stream. If we can determine the direction of the stream producing the dropoff in the 
recoil rate, we have a much stronger basis for associating that stream with a known one. 
We examine the stream direction in the following section. 

B. Stream Direction 

The direction of a stream as well as the speed may be unknown to us. Indeed, a signal 
in a dark matter detector may even be our first indication of some local stream. The 
characteristic time t c is dependent only upon the stream direction, via Eqn. (|17|). so an 
experimental determination of the former can give us indications of the latter. However, t c 
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alone can only give the direction of the stream in the Earth's orbital plane about the Sun, 
insufficient to reconstruct the full stream direction. 

The angle of the stream A from the normal of the orbital plane is encoded in the 
geometrical factor b discussed in Section 111 CI with b = sinA . The parameter b can be 
determined only through modulation effects due to the dependence of v h s (t), Eqn. (fTKj) . 
on the value of b. A small value for b, corresponding to a stream nearly orthogonal to the 
orbital plane, yields only small variations in v h s over the course of a year and modulation 
effects would likewise be small; for a stream perfectly orthogonal to the orbital plane, t> bs 
is constant (b = 0) and there is no modulation. A large value for b, corresponding to a 
stream incident to the Sun along the Earth's orbital plane, yields larger variation of t> b s and 
larger modulation effects. The dependence of the modulation on b manifests itself primarily 
through two effects: the modulation amplitude and the size of the variation in the cutoff 
energy E COjmax — -Eco.min (see Eqn. (J2BJ)), both proportional to b. The modulation amplitude 
is degenerately dependent upon the presumably unknown stream density. The quantity 
-^co.max ~ -^co.min is not, so b can be extracted from a modulation signal (the variation in E co 
is also dependent upon i> , which may be unambiguously determined from the cutoff energy 
as discussed in the previous section). 

From t c and b, the full direction of the stream may then be determined. In Figure 
we show t c as a function of angle in the xy-, xz-, and yz-planes (angle is from the first 
axis indicated in each plane). The characteristic time t c varies most smoothly with angle 
when the plane is near to that of the orbital plane (the xz-plane is the nearest of the three 
shown). The characteristic time t c changes rapidly in the yz-pl&ne (near 30° and 210°) when 
the stream passes near the normal to the orbital plane. The otherwise flat behavior of t c 
in this case demonstrates how t c is independent of the angle A , as the near orthogonality 
of the yz-plane with the orbital plane means rotations in the yz-plane correspond strictly 
to changes in A and not the angle in the orbital plane (except for a 180° phase shift when 
rotating through the normal to the orbital plane). 

The determination of the parameters t c and b from a detector signal is essentially inde- 
pendent of the WIMP mass. Then, even without an understanding of WIMP properties 
(e.g. mass), we can use these two parameters (by converting them to a stream direction) to 
associate a dark matter detector signal with possibly visible structure in the galaxy, such as 
a stellar stream, that is moving in the same directon. Without knowing the WIMP mass or 
the direction of the stream producing the signal in the detector, a dropoff in the recoil rate 
spectrum at a characteristic energy E c cannot be associated with any specific known struc- 
ture (known via other observations), e.g. the Sgr stream, rather than some other unknown 
stream. As described in the previous section, any visible structure could give an indication 
of the speed of the dark matter stream and, therefore, lead to a determination of the WIMP 
mass by the position of the dropoff via Eqn. (|28|). 

C. Determination of Stream Parameters 

As noted by Stiff, Widrow & Frieman [3^ , the modulation signal can be used to determine 
many of the characteristics of such a stream; we note here that such determinations can be 
made mainly from a small energy range about E c and, in the following, outline how to do so. 
The approximate energy of the rapid change in modulation behavior, corresponding to E c , 
yields the stream speed f 0iStr since E c oc t> str . The location of this cutoff energy varies over 
the course of a year; the extent of the variation depends upon the geometrical parameter b 



24 



(the sine of the angle between the stream and the normal to Earth's orbital plane), which 
can thus be extracted. From b and the characteristic date t c (the time of year at which the 
stream moves fastest relative to the Earth; determined from the date of the peak rate), the 
direction of the stream can be determined via the formulas of Section 111 CI The amplitude of 
the modulation can be used to determine the density of the stream. The width of the cutoff 
indicates the velocity dispersion in the stream. While observation of these characteristics 
could potentially be inhibited by limited energy resolutions in dark matter detectors, the 
resolutions available in some current detectors should be sufficiently high as to not pose a 
significant impediment to extracting many of these stream parameters. 

The various parameters, while all potentially extractable from an observed modulation 
signal, will not be equally easy to extract. The characteristic energy E c , characteristic 
time t c , and stream density p str relative to the smooth halo density (assuming the halo 
modulation is observed at energies not near E c ) are the parameters most easily extracted 
from any modulation signal as they have the most dominant effect on such a signal. From 
the first two of these parameters, two components of the stream velocity can be determined. 
The stream dispersion cr str and the geometrical factor b would be more difficult to extract 
as they have less significant effects on the modulation signal; extracting these parameters 
might require much more significant detector exposure and/or a larger detector. The stream 
dispersion will give an indication of the velocity distribution of WIMPS in the stream, but 
the interpretation of this parameter should be limited to an approximation of the distribution 
only as our assumed Maxwellian distribution of Eqn. ((7J) may not be entirely accurate (see 
the discussion at the end of Section II V E|) . Indeed, directional detectors are likely to be 
much more useful in characterizing this velocity distribution. From b, the third component 
of the stream velocity can be determined. Alternatively, if a stream is due to a dwarf galaxy 
being absorbed by the Milky Way, there could be an associated stellar stream that would be 
indicative of the WIMP stream, yielding independent measurements of the stream velocity. 
However, for late infall of dark matter clumps and other possible WIMP stream origins, 
there would be no such stellar stream and WIMP detection alone would be required to 
characterize the stream parameters. 

A stream signal in a dark matter detector may provide a means of determining the WIMP 
mass, as discussed previously. Individual parameters that play a role in this determination 
are discussed here. The parameters E c , t c , and b can be determined from a detector signal 
without knowing the mass of the WIMP. The latter two parameters would give an indication 
of the direction of the dark matter stream producing the signal. Knowing this direction, 
we could associate the dark matter stream with some stellar stream or known structure 
in the galaxy, such as the Sgr stream, and independently determine the speed of the dark 
matter stream (which presumably matches that of any visible, e.g. stellar, components). 
The characteristic energy E c depends only upon the WIMP mass, nuclear target mass, and 
stream speed, so by knowing the velocity of the stream producing a dropoff in the recoil 
rate at E c , the WIMP mass can be derived (the nuclear target mass is well known in most 
detectors). The mass determination from a stream in this manner could be more precise 
than that afforded by the signal from the smooth halo itself. 

Detectors that can determine the direction of WIMP induced nuclear recoils are the 
ultimate goal for characterizing the local WIMP population; however, technology for such a 
detector is still in the development stage. In the meantime, probing the annual modulation 
will also allow us to characterize this WIMP population by, e.g., observing the presence of 
streams in addition to the smooth dark matter halo. 
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VI. SUMMARY 

A dark matter stream presents observational signals unlike that of a smooth background 
distribution (assumed here to be the SHM). The most noticeable difference is the existence 
of a rapid dropoff in the recoil rate around a characteristic energy E c ; this dropoff yields a 
relatively large and non- cosine-like contribution to the annual modulation around E c that 
may be observable even for a stream much less dense than the SHM. Figure Q shows the 
drastic importance the stream can have near E c for the annual modulation. The presence of 
a rapid change in modulation behavior with respect to recoil energy would be an indication of 
a stream or other non-SHM component of the dark matter. Since a stream's effects are mild 
except near E c (and essentially non-existent well above E c ), the presence of a stream should 
not interfere with using the modulation to describe the background distribution, SHM or 
otherwise. Detection of an annual modulation of the cutoff energy together with the annual 
modulation of the overall signal provide a "smoking" gun for WIMP detection. In addition, 
if the WIMP mass is unknown, but a dropoff in the recoil rate can be associated with a 
known stream (known via other observations), the energy at which that dropoff occurs can 
be used to determine the WIMP mass. 
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